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We propose and study various realizations of a Hofstadter-Hubbard model on a cylinder geometry 
with fermionic cold atoms in optical lattices. The cylindrical optical lattice is created by coprop¬ 
agating Laguerre-Gauss beams, i.e. light beams carrying orbital angular momentum. By strong 
focusing of the light beams we create a real space optical lattice in the form of rings, which are 
offset in energy. A second set of Laguerre-Gauss beams then induces a Raman-hopping between 
these rings, imprinting phases corresponding to a synthetic magnetic field (artificial gauge field). 

In addition, by rotating the lattice potential, we achieve a slowly varying flux through the hole of 
the cylinder, which allows us to probe the Hall response of the system as a realization of Laughlin’s 
thought experiment. We study how in the presence of interactions fractional quantum Hall physics 
could be observed in this setup. 


Introduction. - For our understanding of quantum 
many-body systems, considering spherical, cylindrical, 
toroidal, or even more exotic geometries has in many 
situations proven to be of key importance [T]. This is 
because phenomenologically distinct physical properties 
may be revealed by imposing various boundary condi¬ 
tions. Prominent examples along these lines include per¬ 
sistent currents, protected edge states |2], topological 
ground state degeneracies [3], and spectral flow in re¬ 
sponse to fluxes threading the holes of a system |1]. The 
natural question, as to what extent such theoretically 
intriguing constructions can become experimentally vi¬ 
able, can be seen as a challenge in quantum engineer¬ 
ing. A photonic crystal with Mobius strip geometry 
and an artificial flux threading an atomic ring potential 
[6] have been recently realized experimentally. For cold 
atoms in optical lattices d ig with planar geometries, 
remarkable progress has been reported in devising m- 
ng and experimentally probing d uni [HI 112921] gauge 
fields and topological states [221 [23] . It is the purpose 
of the present work to present a microscopic model for 
a cylindrical optical lattice in real space that realizes the 
Hofstadter model [21] of fermionic atoms subject to a 
perpendicular synthetic magnetic held [see Fig. Our 
proposal builds on, and is motivated by advances in gen¬ 
erating light beams carrying orbital angular momentum 
(0AM) |23l 121]. The simulation of artificial magnetic 
fields on non-trivial geometries has also been addressed 
with the experimentally intriguing idea of a synthetic di¬ 
mension |27H32| . i.e. a manifold of internal states coupled 
by Raman lasers that are interpreted as (a small number 
of) lattice sites. In contrast, our work focuses on real 
space cylindrical lattices, which becomes essential once 
interactions are included to realize fractional Quantum 
Hall (FQH) physics [33935] . 

For non-interacting atoms we present a microscopic re¬ 
alization of the fermionic Hofstadter Hamiltonian [21] on 



FIG. 1. (Color online) Lattices with real cylindric (red sites) 
and ring shaped (blue sites) geometry around the focal plane 
of a lens system. Artihcial magnetic fluxes through the surface 
(i&r) and the interior (J?z) [see panels a) and b)j are feasible 
in both scenarios. Panel c) outlines the proposed setup. A 
Gaussian standing wave creating a lattice along z direction 
(black LGoft). Tightly focused lasers carrying orbital angular 
momentum I and —I, respectively (black LGi^p and LG-i^p), 
create the azimuthal lattice potential. Counterpropagating 
lasers with orbital angular momenta h and h (gray LGi^p 
and LGi 2 ,o) create via Raman processes. 

a cylinder 

Ho = Jj e “*■ 

j j 

where J? and denote amplitudes for nearest neighbor 
hopping to lattice site j by a displacement vector z and 
(/)j in axial and azimuthal direction, respectively, oj an¬ 
nihilates a spinless fermion at site j, and the spatially 
dependent phase factors e*'^J affect the neutral atoms 
analogous to the effect of a magnetic field on charged 
particles. Adding interactions, the total Hubbard Hamil¬ 
tonian is given hy H = Hq -\- Hj. For spinless fermionic 
atoms the natural (minimal) interaction is the nearest- 
neighbor interaction Hj = UJ2{ij) where ni = afai 
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is the particle number on site i. 

Below, we address two main questions. First, we detail 
how a cylindrical optical lattice with a synthetic mag¬ 
netic flux <i)r piercing its surface can be created. This is 
achieved by employing tightly focused laser beams that 
carry 0AM [see Fig. Ic) for a schematic], both to create 
the lattice potential and to engineer an artificial magnetic 
held via Raman-assisted tunnelling. In addition, by ro¬ 
tating the cylindrical lattice potential around its axis, 
we are capable of mimicking a slowly varying flux <i>z 
threading the hole of the cylinder in the axial direction 
[see Fig. m (see also Ref. [3]). We argue how the re¬ 
sulting spectral flow hallmarking the quantum Hall effect 
can be experimentally observed. These ingredients allow 
us to implement Laughlin’s original thought experiment 
[3] , explaining the integer quantum Hall effect HJ [33J [37] 
in terms of spectral flow. 

Second, going beyond the single particle picture, we 
investigate numerically the effect of local interactions 
and outline how FQH physics |8flfl3H] can occur in our 
cylindrical model. We compare the Hofstadter-Hubbard 
model on a real space cylinder proposed here to its 
counterpart involving a synthetic dimension [311 |31|. 
Including interactions, synthetic and real dimensions 
are shown to lead to a qualitatively different behav¬ 
ior. The physical reason behind this is that spatial 
locality of many-body interactions and localized wave 
functions giving rise to topological protection occur 
naturally in real samples but may be non-generic or 
require fine tuning in the context of synthetic dimensions. 

Cylindrical optical lattices. - We are interested in cre¬ 
ating optical lattices with annulus or cylinder geometry 
as depicted in Figs[^) and [^). Moreover, building on 
Raman assisted tunneling techniques ndnda and rota¬ 
tion of the lattice potential, respectively, we aim at engi¬ 
neering the two fluxes <I>r $ 2 . The spinless fermionic 
atoms are assumed to occupy a single Zeeman m-state, 
where the quantization axis is defined by homogeneous 
magnetic held B = Bqz. The key ingredients for our pro¬ 
posal are tightly focused light beams carrying 0AM ISi- 
l4(T| . A large focus angle is essential to reach lattice spac- 
ings d on the order of the optical wavelength X = 2TTlk 
and thus sufficiently large energy scales where 

Eh = is the laser recoil energy. Some basic fea¬ 
tures, however, can be understood in terms of a paraxial 
description valid for weakly focused light. There, La- 
guerre Gaussian (LG) laser beams LGi^p with electric 
held ~ describe 

light carrying an 0AM of Ih per photon m- Here is 
the polarization and denote associated Laguerre poly¬ 
nomials and ^ In this framework, an intensity 

modulation I ~ cos^ {l(j)) around the beam axis arises 
naturally when two beams LGpp and LG-pp are inter¬ 
fered |42l - l44] . Together with the radial dependence of 
the intensity, it gives rise to lattice potentials consisting 
of p -|- 1 concentric rings with 21 azimuthal sites each [see 



FIG. 2. (Color online) Panel a): Total intensity pattern of two 
interfering beams: LGrp and LG-7,1 including nonparaxial 
effects, realizing a small two-ring system in the focal plane. 
Panel b): radial cut through the optical potential for an az¬ 
imuthal lattice (interfering LGeo.is and LGeo.-is). Panel c): 
Intersite distances di , d 2 [see a)[ converging to the diffraction 
limit for strong focusing, away from the grey paraxial regime, 
ds refers to the innermost ring in our main setup in [for pa¬ 
rameters see Fig. [^. The upper part of d) shows the natural 
variance of single particle energy levels due to strong focus¬ 
ing. In the lower part it is compensated for by a flattening 
potential (black levels, for details see |60p. and then tilted by 
a magnetic field gradient (blue levels show an exemplary tilt) 
to make two-photon Raman processes between neighboring 
rings resonant. The dotted lines indicate the virtual interme¬ 
diate atomic states involved in the Raman process. 

Fig#)[. Beams with I = 25 and p = 10 have been real¬ 
ized experimentally |45fl49] . The paraxial approximation 
would limit the lattice constant d to values much larger 
than A/2 jSOj. This limitation stands in stark contrast to 
the key requirements of reaching experimentally viable 
energy scales. Here we investigate tightly focused beams 
with 0AM, thus going conceptually beyond the simplified 
picture of the paraxial approximation. The dependence 
of d on the focusing angle 9 |3D] is displayed in Fig. H). 
showing that d « A/2 is achievable in the non-paraxial 
regime. 

To achieve strong focussing, we consider two apla- 
natic lenses with focal length / to focus incident LG 
beams propagating along the z-axis, i.e. the optical axis 
of our setup [see Fig. For an incident LG beam 

with angular momentum I and circular polarization e+ = 
— (e^; -b iey)f^/2, the electric held E/,p close to the focal 
plane is given by ismsn 

E;,, = . (2) 

We note that in this non-paraxial regime the po¬ 
larization can no longer be separated from the 
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spatial mode profile. Moreover, the focused field 
is no longer transverse and the different po¬ 
larization amplitudes are given by (p, z) = 
Eq Jq^'cIO sinO gi {9) V cos 9h„{9)Ji{kp sin0)e*^^ msO ^ 
where we abbreviated h±{9) = 1 ± cos 9 and 

hz{9) = —V^ism9, as well as gi{9) = 
with ^ sin0. The parameters Eq and wq specify 

the intensity and the waist of the beams incident on 
the lens. The numerical aperture of the lens enters via 
NA = sin 9m- The Bessel functions are denoted by J/. 

By superimposing two such laser beams with the same 
propagation direction, but with opposite I, we obtain an 
intensity pattern z), which is invariant under ro¬ 

tation of tt/Z around the optical axis [see Figj^)] in each 
polarization component !„. Importantly, even though 
the azimuthal phase dependence of the electric field is 
not simply [cf. Eq. ([^] like in the paraxial treat¬ 
ment, this symmetry is guaranteed due to the circular 
polarization of the incident light. For details we refer to 
the supplemental material. The optical potential is then 
given by V{r) = -, 7 r}where is 

the atomic polarizability, and we assume = a < 0 be¬ 
low. By adding two Gaussian beams in a standing wave 
configuration we create a lattice along the optical axis. 

In this setting we can realize two different geometries: 

One possibility is an annulus geometry realized by con¬ 
fining atoms to the focal plane [Fig. [^), blue sites in 
Fig-[5)1> consisting of a series of concentric ring lattices 
with d ~ A/2. In the radial direction, the rings are de¬ 
coupled by an energy offset stemming from the radially 
varying laser intensity [see Fig. m- Such an offset can 
be used to engineer the synthetic gauge field by cou¬ 
pling the rings with Raman lasers Hdni. 

Instead we focus here on the cylinder geometry, ob¬ 
tained by restricting atoms to the innermost (energeti¬ 
cally lowest) ring of each disk [see Fig. §)] , and cou¬ 
pling the rings by Raman-assisted tunnelling [a [mill]. 
We note that strong focussing leads to an inhomogeneity 
along the z-direction of the cylindrical lattice potential, 
where the dominant effect stems form the decrease of 
laser intensity away from the focal plane. This leads to 
the trapping potential [see Fig. [^)[ but also to a minor 
z-dependence of the azimuthal hopping m- We observe 
that the latter dependence is even smaller for high radial 
modes p > 0. 

Artificial magnetic field on a cylinder. - To mimic 
the magnetic flux we here adapt the Raman assisted 
tunnelling scheme to our cylindrical optical lattice. This 
scheme consists of two steps. First, the optical lattice 
potential is tilted in one spatial direction, for example 
by a magnetic field gradient, which renders the hopping 
in this direction off-resonant. Second, we restore hop¬ 
ping via a two-photon Raman process. Using LG beams 
(LGq_o and LGi^^, see Fig. [^)) for this Raman process 
allows us to imprint the required phase pattern for the 
flux 

In our setup, strong focusing introduces naturally an 




FIG. 3. (Color online) Laughlin experiment on a cylin¬ 
der. Panel (a): Single particle energy spectrum of Q with 
a threaded flux $2 = Iflyr (black lines) and < 1>2 = 0 (gray 
lines) as a function of the dimensionless azimuthal lattice 
momentum q which is quantized in units of 2-111%/Nm, where 
Nm = 40 is the number of magnetic unit cells. Red cir¬ 
cles show states occupied after flux threading, while blue 
arrows indicate population transfer with respect to the ini¬ 
tial state population (black circles). Panel (b) shows (nj) 
— the average density of the gas after flux $ has been 
threaded (black — $ 2^0 = 0, red — $ 2 ,/ = lArr, blue — 
>I> 2 ,/ = — Idvr). Panel (c): Temperature dependence of the 
visibility u = ((rig —(ug *’'’^))/(’^5 Parameters 

are: maximal azimuthal lattice height F/ = 8.5Eu, axial — 
14 = lEji. The lattice consists of 11 rings each containing 
120 sites (Z = 60,p = 15, h = 20, h = -20). 


optical potential tilted in the axial direction [see Fig. il)l- 
However, this tilt is not a linear staircase, requiring in 
principle separate Raman transition frequencies to cou¬ 
ple neighbouring rings, which would practically limit this 
scheme to a few coupled rings. Gylinders with a large 
number of rings can be achieved by compensating for 
the nonlinear energy offsets by adding an appropriate 
spatially-dependent AG-Stark shift, as familiar from flat¬ 
tening of inhomogeneous optical lattices [S]. By adding 
a magnetic field gradient lllisi in the z direction, we 
achieve a linear tilting as illustrated in Fig. il)- This 
allows us not only to couple the rings by a single pair 
of Raman beams generating a homogeneous flux but 
residual nonlinearities away from the flattened region also 
provide a sharp cutoff for the cylinder in the axial direc¬ 
tion [see Fig. il)l- 

For a Raman pair of lasers carrying 0AM li and I 2 , 
the hopping amplitude between sites j and j -I- z acquires 
the phase dependence 

J.^ = ^ ( 3 ) 

where S M and (/j is the azimuthal angle of the lat¬ 
tice site j. An atom hopping around a plaquette of the 
lattice with 21 azimuthal sites now picks up a nonzero 
phase S(j) = 2stt/1, {s = I 2 — h)- Due to strong focusing, 
different polarization components carry different phase 
dependencies [see Eq. (il- However, these non-paraxial 
effects do not change the phase pattern imprinted 
on the hopping elements. The final tight-binding model 
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is then of the form of the Hofstadter Hamiltonian on a 
cylinder Q . Below we explore its band structure includ¬ 
ing the adiabatic Hall response, which may be probed by 
spectroscopic techniques and time-of flight imaging, and 
the role of interactions. 

Flux threading and Hall response. - The Hall response 
of a cylindrical system can be probed along the lines of 
Laughlin’s famous thought experiment |3], where a small 
electric held in azimuthal direction is imposed on a gas of 
electrons by adiabatically threading a flux in axial direc¬ 
tion through the cylinder. A quantum Hall system such 
as our Hofstadter model 0 then responds with a current 
in axial direction that is generated by spectral flow. The 
integer number of particles transported between the ends 
of the cylinder per threaded flux quantum then equals the 
quantized Hall conductance in units of e^//i. 

Our present setup offers a natural way to realize such 
a scenario. An artificial time dependent flux may be 
implemented by a slowly accelerated rotation with fre¬ 
quency 0(t) of the cylindrical lattice around its symmetry 
axis, which translates the spectrum in quasimomentum 
by Ak{t) = 2mil{t)r^/h close to the focal plane, where r 
denotes the radius of the cylinder. Experimentally, such 
a rotation is readily achieved by a frequency detuning of 
the counter-propagating laser beams generating the lat¬ 
tice potential. At the level of the tight-binding model 
Q, it imprints a complex phase to the azimuthal hop¬ 
ping, i.e., —>■ The small detunings caused 

by the rotation do not affect the other building blocks of 
our proposal. In Fig. 0). we show the Hall response to 
such an axial flux of a system of ultracold fermions with 
magnetic flux = 27r/3 per plaquette [see Eq. 01 at 
zero temperature and 1/3 Ailing of the lattice, i.e., with 
one filled Chern band [see Fig. Depending on the 

orientation of the magnetic flux, atoms are transferred in 
different directions between the edges of the cylinder [see 
Fig. [|3)[. For finite temperatures, the visibility of this 
charge pumping is shown in Fig. §). 

Effect of interactions and fractional filling. - We now 
turn to the role of interactions, focusing on the possibil¬ 
ity of realizing fractional quantum Hall (FQH) physics 
[ 551 - 155 ] on a cylinder. We are particularly interested in 
small cylindrical systems that may be realized in a first 
generation of experiments. Using exact diagonalization, 
we assess which signatures of FQH physics could be ob¬ 
servable. We focus on one-third Ailing of the lowest band 
at flux = 27r/3. In this situation, short-ranged in¬ 
teractions are expected to stabilize FQH states IMUMI, 
more specifically the 1 /3-Laughlin state [341 1551155 ] . In 
our model study we assume a nearest neighbor interac¬ 
tion between the spinless fermions Hi = 
as naturally realized microscopically by magnetic inter¬ 
actions in dipolar Er and Dy fermionic quantum gases 
|56| . We compare our results for the real space cylin¬ 
der to its counterpart involving a synthetic dimension 
m- While the free Hamiltonians are the same and give 
rise to a metallic phase, the interaction in the presence 
of a synthetic dimension is naturally of infinite range 




FIG. 4. (Color online) Azimuthal current jf around the rings 
of the cylinder as a function of the (repulsive) interaction 
strength for the real cylinder a) and the synthetic dimension 
analog b). We consider 6 fermions on a cylindrical optical 
lattice made of 9 rings with 6 sites each (i.e. I = 3,li = 
l,l 2 = —1). Accounting for parity symmetry in 2 -direction, 
only the 5 independent rings are shown. 


|57fl55] along the synthetic dimension, since different lat¬ 
tice sites are not spatially separated in real space, i.e. 

Hsyn = ^ (Sz ^i0-|-zz) ’ 

In the cylinder geometry, a simple probe for a FQH 
state is to look at the interaction-induced emergence of 
chiral edge currents. To this end, we compute the value 
of the azimuthal current it = (ijfa^ - a; -I- h.c.) in the 

^ j+0j J ^ ^ 

groundstate on the different rings 2 as a function of inter¬ 
action strength. As the system is translation invariant in 
the azimuthal direction, the current only depends on 2 . 
Our numerical results are shown in Fig. [^ For the real 
cylinder geometry [see Fig. the groundstate remains 
in the zero momentum sector as expected for an incom¬ 
pressible liquid such as the Laughlin state. At small in¬ 
teraction, the current of the groundstate with local inter¬ 
actions is of the same order on each ring of the cylinder. 
With increasing interaction strength, we observe an in¬ 
creasing dominance of the edge current. While, given the 
small system sizes, the system does not display a fully in¬ 
sulating bulk yet, this may be seen as a first signature 
of FQH physics that could be observed in experiments 
on very small cylinders. Experimentally going beyond 
numerically accessible system sizes, this signature is ex¬ 
pected to become increasingly clear. By contrast, the 
groundstate of the system involving a synthetic dimen¬ 
sion [see Fig. 0)1 has a non-zero momentum for U ^ 0.6 J 
and displays smaller edge currents, decreasing with the 
interaction strength. 

Outlook. - The present proposal can be extended to 
create an optical lattice with torus geometry: this can be 
achieved with two concentrical cylinders which, in our 
setup, are naturally decoupled due to a radial energy 
shift. Coupling of the cylinders at the edges can be re¬ 
stored by photon-assisted tunnelling, effectively sewing 
together the patterns to a seamless torus. 
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SUPPLEMENTRY MATERIAL 

In this Supplementary Material, we present a more de¬ 
tailed and quantitative analysis of the optical potential. 
Moreover, we microscopically determine the relevant in¬ 
teraction parameters for a realistic cylindrical lattice po¬ 
tential. 


I. QUANTITATIVE PROPERTIES OF THE 
POTENTIAL 

The interference pattern defining the lattice structure 
in the azimuthal direction is created by superposing two 
tightly focused copropagating beams that carry opposite 
orbital angular momentum (0AM) li = I, h = —I- For 
the specific example presented in this work we have cho¬ 
sen the Laguerre-Gaussian modes incident on the lensing 
system to have azimuthal and radial indices I = 60 and 
p = 15 respetively. The focusing of the beam is described 
by the ratio of the focal length to the paraxial waist of 
a pre-lens beam, which we choose to be f/wo = 10.35. 
The asymptotic propagation angle of the outermost in¬ 
tensity maximum determines the focus angle of the whole 
beam to the value of 6 = 56.5° [see equations of the main 
text]. 

As discussed in the main text, the total intensity 


I tot ^ \Ei splits into polarization components: 

I = - tt} 0, z), where each of them is of the 

form 

Ia{p, (/>, z) = Qcrip, z) J- bcr{p, z) COS^ (Ifj) J- Ca[p, z)), (4) 

such that Icr(p, <Piz) = la^PiCf) + /I, z). The total inten¬ 
sity /tot = ^ la- satisfies the same periodicity condition 

cr 

and can be expressed as: 

Jtot(p, z) = atotip, z) -I- 6tot(p, z) cos^(Z(/> -I- Ctot{p, z)). 

(5) 

In the ideal paraxial case [1], the polarization of the 
electric field is not affected by the focusing, and in gen¬ 
eral Ua- = Ccr — 0 However, in the nonparaxial regime 
the latter no longer holds. The most important param¬ 
eters of the resulting optical potential are the values of 
Otot J ^tot J Ctot in the proximity of intensity maxima, de¬ 
termining the potential depth, the trapping frequencies, 
and relative shear of the lattice rings. We find that in 
general the parameter c is close to zero, near the focal 
plane. 

To fully control the relevant energy scales, it is im¬ 
perative to have independent control of azimuthal lattice 
depth and radial trapping frequency. To achieve this, 
an intensity imbalance of the lasers LGi^p and LG-i^p 
creating the azimuthal lattice is employed. The resulting 
potential is again of the form 0 , with otot ^ Ii + I-i and 
btot \/hl-i- The independent control of the two co¬ 
efficients, allows to adjust the radial trapping frequency 
(depending both on Otot and btot in 0) independent of 
the azimuthal potential, depending mainly on btot- In this 
work we have assumed that the radial trapping frequency 
of the lattice sites is tuned to make the excitation energy 
of the radial mode almost equal to the excitation in the 
(/) direction (see also Fig.|^. 


A. Band Structure calculation 

When computing the band structure of the total po¬ 
tential we assume that the axial lattice separates the full 
potential into a series of decoupled 2D problems. This 
assumption is valid close to the focal plane, where the 
beam is almost parallel to the optical axis. The resulting 
Schrodinger equation for a single atom in a 2D potential 
reads as 

(-^ A-b U(p,0)) 4'(p,(()) =/;4'(p,(/)) (6) 

In the radial coordinates (energy is measured in units 
of the recoil energy , fc = ^, the length unit 

is defined by setting fc = 1) we obtain: 
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FIG. 5. (Color online) Bloch Band structure for a poten¬ 
tial with ring geometry with ring geometry dehned in terms 
of LGeo.is modes. Here 14 = 7ifij,14 = S.^Er. The en¬ 
ergy is measured relative to energy well depth. The almost- 
degenerate lowest excited Bloch bands are azimuthally ex¬ 
cited (green), radially excited (red). Black dot marks the 
well-defined lowest Bloch band, while blue dots denote second 
azimuthally-excited bands. Soft radial modes spawn a whole 
ladder of states characterized by azimuthal excitation. Exact 
formal identification of "radial" and "azimuthal" excitations 
is not possible due to weak p — <j> coupling. 
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FIG. 6. (Color online) Panel a) shows the axial hopping am¬ 
plitudes jt for different rings of the lattice potential. In 
this setup other hopping directions are prohibited by en¬ 
ergy conservation. Panel b) shows the two-body interac¬ 
tion integrals L(j in the respective sites. In both panels the 
dashed red frame marks lattice sites that are part of the 
cylinder considered in the main text (sites for z < 0 are 
not shown due to symmetry). The potential parameters are: 
f/wo = 10.35, = 8.5,14 = 7Er. 


We want to find the Bloch states '0,, 


(P,'/') ='«5(p,0)e*-‘^, (8) 

where u is a j periodic function in the 0 direction. Ex¬ 
panding i/jg, in a Fourier series in the 0 variable we get 


Uq{p, 0 ) 


1 

V2l 




m——oo 


(9) 


where q is the dimensionless quasi-angular momentum. 
The 2D Schrodinger equation Q becomes then an eigen- 
problem expressed as a system of coupled differential 
equations in ID real space. We solve these equations 
using sparse matrix diagonalization of a discretized real 
space equation. 

The formal singularity of the eigenproblem (|^ at 
p = 0 is resolved by imposing proper boundary condi¬ 
tions. Indeed, all components of Eq. ([^ corresponding to 
a nonzero angular phase dependence must vanish at the 
origin, while the only component with no angular phase 
dependence, do{p,0) must satisfy: 0pdo(Pj 0)|p^o “ 9. 

The numerical solution of the eigenproblem given by 
Eq.j^gives rise to the Bloch spectrum shown in Fig.[^ It 
features a lowest Bloch band separated from all excited 
states. As the problem Eq. is weakly-nonseparable, 
the overall structure of the excited Bloch states in the 
spectrum is well described by elementary radial and az¬ 
imuthal excitations, coupled by very small avoided cross¬ 
ings. 


B. Tight-binding and Hubbard parameters 


The purpose of this section is twofold. First, we mi¬ 
croscopically compute the hopping integrals for the 
cylindrical optical lattice. Second, while the focus of the 
discussion in the main text is on spinless fermions, we 
here present the on-site Hubbard interaction parameters 
of the cylindrical lattice potential. These parameters are 
relevant if the cylindrical optical lattice is loaded with 
bosons which is experimentally equally conceivable. 


From the numerical solution of Eq. we can extract 
a set of Wannier functions in a standard way |2].The 
Hubbard parameters describing the contact interaction 
Uj = ^^^/dV|wj(7=0|'‘ and the hopping integrals jf 
are determined from band computation and construction 
of the Wannier functions. 


The t/j, parameters are constant within each ring 
formed by lattice sites due to rotational symmetry. 


In Fig. we show curves denoting the U, J parameters 
for different potential rings. Different curves show val¬ 
ues for the innermost, second-innermost etc. rings. We 
note here that the azimuthal lattice depth in the inner¬ 
most ring is the deepest one, which results in the lowest 
hopping rates shown in Fig. For larger ring numbers 
the increase of radius and the corresponding distance be¬ 
tween the lattice sites eventually leads to a decrease of 
the hopping rate. 



































FIG. 7. (Color online) Panel a) shows effect of rotation on 
energy levels of band computation. In addition to the shift 
of the dispersion relation in the dimensionless quasi angular- 
momentum q ^ q + 2mQ,r^ Ih it introduces also an irrelevant 
offset of the energy /S.E. Panel b) shows dependence of flux 
^z,n threaded by imprinting the phase (j>z on azimuthal hop¬ 
pings on the ring number n in cylinder configuration. Colors 
Black/Red indicate fluxes $z,o, ^z ,5 as a function of II. Rings 
indicated blue are beyond the range of the cylinder considered 
in the main text. 


II. MICROSCOPIC MODELLING OF 
ARTIFICIAL MAGNETIC FIELDS 

A. Rotation of the optical lattice 

The flux through the cylinder is implemented by 
rotation of the lattice in the azimuthal direction. A slight 
difference of frequencies of lasers creating the azimuthal 
optical lattice oji — 0 J 2 = causes its rotation around 
the optical axis with angular velocity [T]. The time 
dependent transformation 

( 10 ) 

leads in the co-rotating frame to a modified Hamiltonian 
as defined by the Schrodinger equation 

ihOti) = (^H + Ip. ( 11 ) 

At the level of the tight-binding description, the ex¬ 
tra term just amounts to a shift in the angular quasi¬ 
momentum. The corresponding shift of the Bloch bands 
is shown in Fig. Formally, this shift is reflected in 
a phase imprinted on the azimuthal hopping amplitudes 
Jf —>■ The phases (pz,n are slightly different for 

different potential disks. For the specific set of param¬ 
eters presented in this work the dependence of (pz,n on 
disk number n is presented in Fig The rotation cor¬ 
responds to threading an artificial magnetic held in 2 : di¬ 
rection, and therefore the flux threaded through different 
rings depends on their radii. 

For the simulation of the Laughlin experiment [see 
Fig. 3 in main text], we consider in this work a cylinder 
with an axial length of 11 sites around the focal plane 
(with site no. 6 being in the focal plane). In this case, if 


the phase imprinted on a lattice ring in the focal plane is 
‘J’z.o = IdTT, the phase imprinted on the outermost rings 
(5th nearest neighbours) is $^^5 « 1 . 02 <i )2 q- Such a small 
variation of the probe flux does not have qualitative ef¬ 
fects on the Laughlin experiment. 


B. Raman-induced hoppings 


A Raman-assisted tunnelling scheme based on a sin¬ 
gle pair of lasers, requires a linear tilt of the optical lat¬ 
tice in axial direction. In our setup this requires two 
steps: 1 ) flattening of the natural single particle offsets 
Aj, 2 ) application of a linear gradient by e.g. a magnetic 
field. 

Regarding 1), the natural strong variation of the single¬ 
particle on-site energy offsets Aj can be flattened using 
external Gaussian beams, as employed to flatten inho¬ 
mogeneous optical lattices |5]. The necessary correcting 
potential should ideally depend only on the z coor¬ 
dinate, and should vary over the distance of a few wave¬ 
lengths. 

These properties are realized e.g. by a potential cre¬ 
ated by Gaussian beams propagating in the perpendicu¬ 
lar direction to the optical axis, focussed at the the same 
point as the lasers creating the lattice. In general, such 
a beam, in the paraxial description, propagating along 
the X axis, gives rise to the optical potential: 




W 

Wy{x)Wz{x) 


exp 





with Wi{x) = woy/l + {x/xr)'^,xr = WqTt/\. To mini¬ 
mize the residual variance of Ej +Vc on, one can focus it 
by a cylindric lens in the x — y plane, i.e. —)■ 00. 

An alternative energy correction scheme could be 
based on strongly focusing correction fields (again of 
Laguerre-Gaussian type) using the same optical setup as 
for the lattice lasers. To achieve maximal dependence 
with the z coordinate, a beam with p = 0 should be 
used. This approach retains perfect rotational invariance 
of the correction held on the innermost ring of the lattice. 

In what follows we have simulated an exemplary Gaus¬ 
sian correction scheme (the same parameters as those 
used in Fig. 3. in the main text) [Eq. @1 • The opti¬ 
mization of parameters of correction lasers is performed 
as follows. First, by a tight-binding computation, we de¬ 
termine the single-particle energy levels Aj of the cylin¬ 
dric potential as a function of ring number. As a cor¬ 
rection potential, we use two Gaussian potentials, one 
with paraxial waist wi, the other with 102 - Then we de¬ 
fine the following auxiliary function with free parameters 

a,/3,uii, W 2 : 


F{a,/3,wo,wi) = 

i^R-1 

= ^ \cz+i-Czf-^\ER+,-En-6\, (13) 

i—L 

Ci = El -aVc{wi, i) - I3Vc{w 2, i) (14) 
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FIG. 8. (Color online) Time-of-fUght images of a lowest 
band Bloch function in a cylindrical optical potential, released 
abruptly from the innermost ring of a potential in the focal 
plane. Here 14 = 7Er,V^ = S.SEr, f/wo = 10.35. Pan¬ 
els show images of different quasi-angular momenta q. Subse¬ 
quent panels show: a) g = 15, b) g = 50. Remaining panels 
show g = 50 as well in the case when the optical potential is 
switched off adiabatically, to allow for adiabatic band map¬ 
ping HI- Panel c) deals with the case when adiabatic band 
mapping is performed by ramping down just one of laser cre¬ 
ating the lattice, thus reducing 14 >, from S.bEu to a final value 
0.5Eii, while in panel d) intensity of both lattice lasers is ho¬ 
mogeneously decreased to the point when V), = 0.5Eji. 


and Lagrange multiplier 7 . Minimization of F ensures 
that Ei are flat between sites L = — 5 and R = 5, and 
that the energy S barriers deflning boundary condition 
are high, allowing a hard wall description. Here we fix 
(5 at (5 « 5.5J^. The energies of single-particle states of 
particular rings are then almost equal Ei « E+AEi with 
I AS, I < 0.08 JL 

Microscopic Raman-assisted hopping parameters. 
With the application of a linear tilt of the onsite energies 
of A per lattice constant the bare tunnelling elements in 
axial direction are off-resonant, A. Tunnelling 

can then be restored via a two photon Raman process. 
This allows also to introduce the flux As discussed 
in the main text, to create <i>rj requires a pair of counter- 
propagating Laguerre-Gaussian beams. The Raman- 
induced tunnelling amplitude between two s-band states 

is given by J? « Jf’^'^'^Ji{2^rizr]^) |3], where denotes 
the height of the time-dependent potential created by the 
Raman lasers. The quantities rjz = (wj| exp( 2 iA; 2 ;)|xi;j+z) 
and rpf, = (wj | exp(|^(/))|z(;j_|_^) are defined via integrals 
involving the Wannier functions wy For our setup, 77 ^ 


and Qj are j-dependent due to focussing-induced varia¬ 
tion of the azimuthal lattice depth. However, the vari¬ 
ance of 770 is small, and the Raman lasers, operated far 
from the diffraction limit, exhibit only a minor intensity 
variation over the innermost rings. In conclusion, our 
microscopic calculations show that the resulting varia¬ 
tion in the Raman-induced hopping is less than 1% over 
the whole cylinder and is thus much smaller than the 
variation in the natural azimuthal hopping. 

III. TIME OF FLIGHT IMAGES 

Time-of-flight (TOF) imaging is a standard experimen¬ 
tal method to probe a single particle Green’s function 
G{r,r') = {'4){r,t)'iJ)\F,t)) (see |4]). After the atomic 
cloud is abruptly released from the confinement and al¬ 
lowed to expand freely, an absorption imaging is per¬ 
formed along the optical axis z which measures the den¬ 
sity of the gas after time ixoF- The measured den¬ 
sity is the full 3D density of the cloud, i.e., 7 r(x) = 

l^(k)P^(k = integrated along 2 : di¬ 

rection, where w is the Wannier state of the occupied 
Bloch band. 

Fig. [8] shows the absorption image for a single Bloch 
state calculated using the microscopic description of the 
potential modelled in this article (with I 4 = ^En, = 
8.5i?fl, f /wq = 10.35), and a fixed quasi-angular momen¬ 
tum q. It has the structure of a ring with a radius propor¬ 
tional to |g|. This allows to measure the distribution of 
the absolute value of the angular quasi-momentum. Ro¬ 
tating the lattice as described in |H A| allows to shift the 
quasimomentum distribution and measure it as a func¬ 
tion of |g —go| with the experimentally controllable offset 
Qo ■ 

An alternative approach is to moderate the speed with 
which the lattice is switched off before the free expansion 
phase (!]. This enables the Bloch states to adiabatically 
follow the decreasing potential height. This converts the 
initial quasi-angular momentum for the deep lattice to 
the ordinary angular momentum distribution. This al¬ 
lows to reduce higher order peaks for large angular quasi¬ 
momenta q as seen in Fig. |^). 

We consider two cases of adiabatic ramping down of 
the potential. In the first case [see Fig. |^)] only one of 
the lasers creating the azimuthal lattice is ramped down 
(and the other is slightly ramped up to keep the radial 
curvature of sites constant). In terms of numbers the 
initial azimuthal lattice depth = S.BEn is reduced to 
— 0.5. In the second case [see Fig. I^)[ both lasers 
creating the azimuthal lattice are rampea down until the 
depth of the azimuthal lattice reaches a final value of 
= 0.5. Afterwards the lattice is switched off abruptly 
(close to 14 « 0 , adiabaticity is not possible). The Lasers 
creating the axial lattice are also switched off in both 
scenarios. With this technique the best discrimination 
between different jgj can be obtained in the TOF image. 
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